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I analyze non-local symmetries of finite-size Euclidean 3D lattice Chern-Simons models in the 
presence of an external magnetic field and non-zero average current. It is shown that under very 
general assumptions the particle- vortex duality interchanges the total Euclidean magnetic flux $/27r 
and the total current I in a given direction, while the flux attachment transformation increases 
the flux in a given direction by the corresponding component of the current, ^/2n -^ $/27r + 21, 
independently of the disorder. In the language of 2-f 1 dimensional models, appropriate for describing 
quantum Hall systems, these transformations are equivalent to the symmetries of the phase diagram 
known as the Correspondence Laws, and the non-linear current -voltage mapping between mutually 
dual points, recently observed near the quantum Hall liquid-insulator transitions. 
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INTRODUCTION 

Beginning with the |-analysis of the Ising model by 
Kramers and Wannieild in early forties, the particle- 
vortex duality found extensive use in physics. However, 
duality had been always considered as a kind of "the- 
oretical" symmetry, whose manifestations are limited to 
the critical region and generally are very hard to jobserve. 
This is one of the reasons why the experimentg3'BQ ob- 
serving a precise duality between the adjacent quantum 
Hall phases came as a total surraise. The agreement 
with the apparently simple theorycl combining the dual- 
ity transformation and the random phase approximation 
(RPA) went far beyond the original expectations of the 
authors. In fact, experimentalljcl the duality was shown 
to relate the entire set of non-linear I-V curves in a sur- 
prisingly wide region of temperatures, gate voltages and 
magnetic fields, indicating that the dynamics on the two 
sides of the transitions must be very similar. This is 
in spite of the fact that theoretically, according to the 
very same duality mapping, the elementary particle ex- 
citations at the two sides of the transition, fermions and 
fractionally charged vortex-like quasiparticles, are very 
different! Even more surprising are the experimental in- 
dications of duality betweeHiihe two phases in the metal- 
insulator phase transitionsBQ, apparently observed in sil- 
icon MOSFETs in the absence of a magnetic field, where 
the conventional theory of localization predicted no phase 
transition at all. 

The fractionallypi charged quasiparticles were intro- 
duced by LaughlinH as excitations on top of the highly- 
correlated state now known as the Laughlin liquid; these 
objects were used to explain the observed fractional val- 
ues of the Hall conductance. Unfortunately, the wave- 
function description does not easily reveal the order pa- 
rameter associated with different phases and it could 
not be used to study the phase transitions between 
the plateaus. The interpretation of the peculiar off- 
diagonal long-range order inherent to Laughlin's wave- 
function as the condensation of the composite particles 



was proposeda by GLcYin and MacDonald. These ideas 
were later developedtd by Zhang, Hansson and Kivel- 
son, who derived the Chern-Simons-Landau-Ginzburg 
(CSLG) theoryM directly from the microscopic Hamilto- 
nian. The symmetries of the ground-state wavefunction 
now became apparent as the (approximate) non-trivial 
symmetries of the model itself. These symmetries impose 
many important constraints on the allowed phase transi- 
tions and were usedEl by Kivclson, Lee and Zhang (KLZ) 
to construct the Global Phase Diagram of the quantum 
Hall system. 

Despite all the success of the CSLG,|theai:y,|aBd| its ex- 
cellent experimental confirmationEj'llaEJllaoHlj'Q, some 
questions still remain unanswered. The most important 
of these questions is why does so simple a model work 
so well. The original KLZ calculations were done only 
within the RPA approximation, and the universality of 
the scalar field correlators was conjectured rather than 
derived. Indeed, the quasiparticles have different charges 
and even their fractional statistics differ, and one would 
think that the global universality in this pictm-e is not 
likely to happen. For example, calculationgi3'E3, pertur- 
bative in the Chern-Simons coupling, indicate that this 
additional gauge coupling is at least a marginally rele- 
vant pcrtupbation to the model of free bosons. Similar 
calculationli3, performed without the 1/iV expansion sup- 
pressing the effect of gauge degrees of freedom, demon- 
strated that the Chern-Simons gauge interaction is actu- 
ally a relevant perturbation to the scalar field fixed point; 
this interaction leads to a first-order phase transition in 
the regime where perturbation theory applies. 

A formal symmetry-based approach to the proMcHi of 
the quantum Hall phase transitions was initiatfidE2[Eil by 
Lutken and Ross, and later further developedcjS to ap- 
ply for the beta-function flows. The idea was to find 
all the consequences of the observed SL(2,Z) symmetry 
group without specifying how this symmetry is achieved. 
In the spirit of the j-twa-parameter scaling model of the 
quantum Hall cffcctcJ'Ej, it was postulated that the par- 
tition function and the scaling beta-function for the com- 
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plexified conductance S = H^y + iT,xx of the system, are 
invariant under the elements of S'L(2,Z) and can depend 
on E only as "quasi-holomorphic" functions. This as- 
sumption, together with the known behavior in the lim- 
its of weak and strong coupling, severely restricts the 
phase diagram and the possible functional form of beta- 
function. Unfortunately, so far nobody has managed to 
complete this program and find the critical exponents in 
the vicinity of the critical point relying on the symmetry 
considerations, alone. 

As an alternative to the conventional perturbative ap- 
proach, S.-C. Zhang and the author attemptedc3 to un- 
derstand the universality of the quantum Hall transi- 
tions using the lattice x-y model coupled with the Chern- 
Simons gauge field, which comes naturally as the strong- 
coupling regime of the CSLG model. They constructed 
the exact duality and fiux attachment transformations 
on the lattice and proved that the Chcrn- Simons field is 
the most relevant interaction in this problem and, mod- 
ulo the irrelevant higher-order in momenta terms, the 
phase transitions in different Chern-Simons models with 
the coefficients related by these transformations mus4 be 
in the same universality class. Fradkin and KivelsonE3, in 
addition to the model similar to that of Ref. I2I, consid- 
ered a model where the gauge-field polarization operator 
acquires an additional T-symmetric, linear in momenta 
part. The corresponding dimensionless coupling con- 
stant, combined with the usual Chcrn-Simons coefficient, 
forms a complex-valued parameter z which transforms 
according to a subgroup of the modular group SL(2,Z). 
In particular, there are infinitely many fixed points that 
remain invariant under certain elements of this group. It 
was proposed that corresponding values of the parameter 
z must coincide with the critical points of the quantum 
Hall system. Once again, symmetry requires the conduc- 
tivity in these critical points to have certain (generally 
different for different fixed points) universal values. 

The major advantage of Chern-Simons x-y models is 
that they allow to treat the gauge fields, apparently the 
most relevant part of the interaction, non-perturbatively, 
considering all other interactions as less relevant. On 
the other hand, the lattice models considered so far did 
not have precisely the sanie-|Symmetries as the contin- 
uum CSLG model derivedllj directly from the micro- 
scopic Hajailtonian of electrons. Specifically, the lattice 
modelsEj^cJ have the exact rotational symmetry (or rela- 
tivistic one after the Wick rotation) while the scalar field 
in Refs. |lO|, ^ is non-relativistic. In addition, none of the 
lattice models was constructed carefully enough to allow 
the finite-size scaling analysis in the presence of external 
magnetic field, current or charge density. 

The goal of this work is to construct a lattice Chern- 
Simons model as a regularization of the continuum CSLG 
model in the limit of strong short-range electron-electron 
repulsion, at the same time keeping track of the long- 
distance properties of the model. I show that the required 
model is exactly the Chern-Simons x-y model, extended 
to include a non-zero flux of external magnetic field and 



a non-zero average current density. Even though this 
model is written in rotationally (relativistically) covari- 
ant form, the Lorentz symmetry is broken by the external 
current (the last component of which corresponds to the 
charge density) and the external magnetic field, so that 
both by construction and by its symmetry the model rep- 
resents a quantum Hall system. It is shown that under 
very general assumptions, in addition to transformations 
of the Chern-Simons couplingn^, the duality interchanges 
the total magnetic flux $/27r and the total current I in 
a given direction, while the flux attachment transforma- 
tion increases the total magnetic flux in a given direction 
by the corresponding component of the total electrical 
current, */27r -^ ^/2tt -|- 21. 

The derived transformations, after being reinterpreted 
in terms of the 2-1-1 dimensional systems, are equiva- 
lent to the symmetries of the phase diagram known as 
Correspondence Lawa3, as well as the nonlinear current- 
voltage mapping between mutually dual points, recently 
observeda near the quantum Hall liquid-insulator tran- 
sitions. Even though the exact symmetries of the phase 
diagram can be derived within othpt approaches, like, for 
example, the 2D GET formulationca, only within the CS 
x-y model one can understand the non-linear transport 
relationships as well. Unlike the previously studied con- 
tinuum models where the problems of renormalization 
arc unavoidable, here the current-voltage mapping is ex- 
act and is an immediate consequence of the duality re- 
gardless of the specific form of the gauge and scalar-field 
coupling. Moreover, these symmetries hold in the pres- 
ence of arbitrary potential and magnetic disorder, which 
remain unchanged by either the duality or the flux at- 
tachment (up to a short-radius formfactor coming from 
the different size of the quasiparticlcs in different rep- 
resentations). This implies that the derived results will 
hold even if the disorder, marginal by the power counting, 
turns out to be a relevant perturbation. 

The Section m introduces the lattice Chern-Simons x-y 
model in an external gauge field as the lattice regulariza- 
tion of the phase part of the CSLG model after integra- 
tion over the density fluctuations. The infrared behav- 
ior of this model is controlled by defining it on a torus 
using the magnetic translation group. The duality and 
fiux attachment transformations of the lattice model are 
discussed in Sections || and III respectively. The appli- 
cation of derived results to the real-time quantum Hall 
systems, including the derivation of the Correspondence 
Lawaj and the current -voltage duality, are presented in 
Sec. 



IV 



Before ending this Introduction, I would like to remark 
that there have been a great number of publications on 
the problem of phase transitions and the duality in the 
quantum Hall system, and it would be impossible to list 
even all the approaches here. An excellent, summary of 
plausible interpretations of the experimenta and the key 
difficulties associated with these interpretations was re- 
cently givenEj by Shimshoni, Sondhi and Shahar. 
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I. DEFINITION OF THE MODEL 

A. Lattice and continuum Chern-Simons models 

The Euclidean partition function of the Chern-Simons 
x-y model with an external field can be written as 



Z ^Y[ I daddt 



-iSe{0,si+A)-iSa{a) 



where the gauge part 
and scalar part 



d^k 



a_k/C(k)ak, 



i5e = -— ^ cos(V6' - a - A) 



(1) 



(2) 



(3) 



of the action depend on the phases 9 defined at every 
vertex of the three-dimensional cubic lattice, while the 
Chern-Simons gauge field a and external gauge field A 
live on the links. As usual for lattice gauge theories, the 
value of A on the link can be interpreted as the integral 
of the physical gauge field along this link 



^n,n+/i 



A^(x)dx'^; 



for notational convenience, in this paper the contour 
sums of gauge quantities are denoted as the continuum 
integrals along the appropriate contours. Similarly, the 
flux of the magnetic field through some area, a sum over 
plaquettes in the discrete lattice model, will be denoted 
as the continuum surface integral of the field over this 
area. These notations were chosen for the sake of read- 
ability of the paper. The Reader interested in the details 
of the lattice definition of the Chern-Simons term is di- 
rected to Ref. ^. 

The dynamics of the field a is determined by the trans- 
verse gauge kernel /C(k). The most general rotationally- 
invariant form of this kernel can be written as the sum 
of T-antisymmetric and T-symmetric parts 



/C^'^(k) = /(fc)e^ 



9ik) 



p2g^.. _ pp 



(4) 



depending on the lattice momentum P^ = 2sin(A:p/2). 
ft is the small-momentum behavior of the functions / 
and g that determines the properties of the model, as 
has been shown in Refs. |2^, pTl Specifically, the usual 
Chern-Simons model can be defined by the condition 



lini f(k) = — 



f 



2-n a 



and 



lim kq(k) 



0. 



(5) 



The generalization oij-ihis standard form considered by 
Fradkin and KivelsoncJ, has both limits finite. 



lim /(fc) = ^ 



and 



lim kq(k) 



go 

27r' 



(6) 



Unlike the usual Chern-Simons model, the model with 
such a kewiel has self-dual fixed points; they were 
associatedE3 with the fixed points of the quantum Hall 
phase transitions. 

As defined, the partition function (Q) is symmetric with 
respect to simultaneous change of the sign of all compo- 
nents of the fields V6', a and A; since the phase gradient 
is associated with the current, both the external magnetic 
field and the current are reversed. 



J, a 



a. 



A^ -A. 



(7) 



This trivial symmetry preserves the handedness of the co- 
ordinate system, and therefore does not change the sign 
of the pseudoscalar a. One can also change the handed- 
ness of the coordinate system, simultaneously reversing 
the sign of a and all pseudovectors like the magnetic 
field. In terms of the global flux of the total magnetic 
field # — / V X (A + a) ds, this transformation can be 
written as 



$ 



*', 



(8) 



while polar vectors, (for example, the particle current 
density j), remain unchanged. 

At first glance, the lattice model (|l|) seems to be only 
vaguely related to the quantum Hall problem. As de- 
fined, this model treats aU|[GOordinates evenly, unlike the 
continuum CSLG modeltSllil with the action 



D — Oa ^" Oir 



did X La+ \ did x £, 



(9) 



which has a similar gauge part La = e'"'''a^9i/ap/(47ra), 
but an explicitly non-relativistic scalar part 



C. 



ip^ {idt 
I 
2mp 



at - At) if 
iV + a + A] if 



(10) 



' ^( t -^2 



defined in terms of the wavefunction ip = ^fp exp id oi the 
bosonic condensate, where 9 and p are the corresponding 
microscopic phase and density. 

Note, however, that even in the absence of the gauge 
fields, the lattice model (ll) represents only the phase sec- 
tor of the continuum model (|9|) ; any such model should 
have a linear-spectrum mode as long as there is a fi- 
nite phase stiffness generating the propagation velocity 
V. Formally, the relativistically-symmetric form of the 
phase action can be obtainedEJ by integrating away the 
fiuctuations of the density p — p^p near its average value 
p, which leads to the Lagrangian 






^{dt0 + at+ Atf - (V0 + a + a)' 



(11) 



dependent on the phase only. Here k — 'p/m,, and v"^ = 
'pV/me are the effective stiffness and velocity of the phase 
excitations calculated with frozen gauge fields; certainly. 



"Global symmetries in quantum Hall systems: lattice description" by Leonid Pryadko 



these values should be understood as bare values of the 
parameters. 

Even after the density integration, the "relativistic" 
symmetry of the model ( [ll|) is not exact: it is broken 
by the external gauge field A which rcficcts the pres- 
ence of a uniform magnetic field and both scalar and 
magnetic disorder. After the Wick rotation t -^ ir, 
At — > —iAz, the normal component of the magnetic field 



B = d^Ay 



^y-^x 



becomes the last component of the 



Euclidean magnetic field B^ = B, while the components 
of the lateral electric field 



E 



-VAt - dtA 



correspond to the remaining components B^ = ie^^Ei] 
clearly, the analytic continuation from the imaginary val- 
ues of external fields is needed to relate the results of 
the two approaches. Similarly, physical charge density 
p and the current density J determine the three compo- 
nents of the Euclidean current density j'^ whose average 
values can be obtained from the partition function dl^) 
by differentiating with respect to the components of the 
Euclidean gauge field, 



ISlogZ 
i SA,, 



(12) 



This implies the correspondence j^ — > (p, if); the spatial 
components of the Euclidean current acquire imaginary 
value and, once again, the analytic continuation must 
accompany the Wick rotation. 



B. Magnetic translations and periodic boundary 
conditions 

By going from the continuum Lagrangian ( |ll| ) to the 
lattice action (0) we regularized all possible ultra-violet 
divergences. To finish the definition of the model, it is 
now necessary to regularize its infrared behavior. This 
can be done by defining the model on a closed surface 
like a torus. Let us imagine the finite-size system of in- 
terest as a unit cell of an infinite lattice, where all gauge- 
invariant quantities are perfectly periodic. On the other 
hand, a gauge field cannot obey the periodic boundary 
conditions in any two directions if it has a non-zero mag- 
netic fiux in the perpendicular direction. The exact pe- 
riodicity in such a system can be established using the 
operators of magnetic translation, defined as the combi- 
nation of the ordinary translation and the special gauge 
transformation restoring the original form of the Hamil- 
tonian. 

As an example, consider a simplified version of lattice 
model (H) where the Chern-Simons field has been frozen. 



iS„ 



i^cos(V0-A) 



(13) 



which is just the three-dimensional x-y model minimally 
coupled to an external magnetic field. Even though the 
vector potential changes under the gauge transformation 



' + /, A ^ A + V/, 



(14) 



the magnetic field B = VxA and all physical observables 
remain invariant. 

Similarly, a displacement by a period L^ in a given 
direction fi — x, y or z changes the vector potential by 
the increment defined as A A(r) = A(r + L^) — A(r). 
Since the gauge-invariant magnetic field VxA is periodic 
by assumption, its increment A vanishes, A (VxA) = 
VxA„A " - - . . _ . 



0, so that the increment A of the vector 



potential A can be expressed as a gradient. 



A A(r) = VF^(r) 



(15) 



The newly-defined scalar function generates the gauge 
transformation dlj) with / = — i^^, which exactly re- 
stores the original form of the model. Therefore, even 
though the exact form of the model (O) is changed by 
both the translations and gauge transformations, their 
combination preserves each term of this model. Such 
combination of a translation and a gauge transforma- 
tion define the "magnetic translation" r(L^). As in the 
case of uniform magnetic field, the magnetic translations 
r(L^) and r(L^) commute with each other only if the 
magnetic fiux through the area L^ x L,y is equal to an 
integer number of flux quanta (fiQ — hc/e. 

Let us utilize the magnetic translation symmetry to 
specify the periodic boundary conditions for our lattice 
system in the region < r^ < L^. The angle 9 in the point 
(x^y^hz) outside the boundaries of this region must be 
replaced by the appropriately corrected angle in {x, y, 0) 
using the symmetry with respect to magnetic translations 



,y,L^) ~ Fz{x,y,Lz) = 0{x,y,O)] 



(16) 



similar patching conditions must be defined for the two 
remaining "outer" surfaces of the system. The ex- 
pression (nw contains only the part of the function 
Fz{x,y,'Lz)a,t the "top" face z = L^ of the this re- 
gion of interest; similarly only the values F^lh^, y, z) and 
Fy{x, hy, z) of two other patching functions at their cor- 
responding "upper" faces are relevant. Each of these in- 
plane values obeys a simpler two-dimensional equation 



A^Aj 



Vif, 



(17) 



instead of Eq. (|T^) , where the index "_L" denotes the part 
of the vector orthogonal to the axis "/i" . In the following 
it is assumed that the patching functions F^ are functions 
of only two variables in the corresponding plane r^ = L^, 
/i = a;, y, or z. 

Generally, under the gauge transformation (|lj) the 
fields A at the opposite faces of the system transform 
differently, so that the values of the increments A Aj^ 
change. To preserve the condition (|l^), we need to adjust 
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the patching functions, Ff^ — > F^ + A /. This transfor- 
mation rule, along with the original definition (|l7|), can 
be satisfied by defining 



replacing the summation over m by the summation over 
new independent integer-valued variables j also defined 
at every link of the lattice. 



F^ir^) 



A(r) • dr + Const, 



(18) 



BCDA 



where the integration is performed along the three-part 
contour starting in the point B = r^ + L and ending 
in A = r^ as shown in FIG. |l|: two congruent lines BC 
and DA along the opposite faces of the system, connected 
by a fixed contour CD entering both faces in equivalent 
points. An arbitrary gauge- invariant constant introduced 
in Eq. ( |l8| ) accounts for the freedom of uniformly "twist- 
ing" the boundary conditions for 0; it is also the freedom 
of choosing the configuration of the fixed segment CD of 
the "back flow" contour BCDA. 
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FIG. 2. Solid lines denote the definition region of the 2D 
Villain x-y model, while the filled circles show the correspond- 
ing vertices of the dual lattice. Region outside the thin solid 
contour represent the images of the points after translation (s). 
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FIG. 1. The coupling of the gauge field with the current 
line AB is made gauge-invariant by introducing the surface 
patching term (llq), equivalent to integration along the con- 
tour BCDA. The line BC is the image of AD after translation 
by the period L^ . 



Normally, the term with the gradient of 6 in Eq. ( poj ) 
can be integrated by parts, ~ij ■ V6 -^ iOV • j, so that 
the subsequent integration in 9 renders the conservation 
law 



Vj = 0, 



(21) 



which, given the form of its coupling to the gauge field, 
allows to identify j as the conserved current. In the finite 
system, however, the phases 6 in equivalent points across 
the boundary of the system are related by the patching 
equations (ttq), and the integration by parts generates 
additional boundary terms 



C. 2D example 

Before advancing any further the analysis of the three- 
dimensional problem, let us analyze a simpler two- 
dimensional Villain x-y model 



^'5v = ^E(^^-A 



27rm) 



(19) 



in external magnetic field. In expression (n9) the angles 6 
are defined in the vertices of two-dimensional square lat- 
tice, the vector potential of the external magnetic field 
is defined at the links, as well as the additional integer- 
valued field m which restores the periodicity broken after 
changing the cosines of the original x-y model for parabo- 
las in the Villain approximation. 

At every link the summation in appropriate compo- 
nent of m is independent of anything else and can be 
transformed with the Poisson summation formula. 



— (V6I - A - 2Timf 



^-zj.(V0-A). (20) 



ZOsi 



-* / ^ r-F^{y)dy~i [ ' f-Fyix)dx, (22) 











where the patching functions Fi are related to the incre- 
ments of the gauge field across the system. 



A^Ay = VyF^{y), A A^ = V^Fy{x) 



(23) 



As before, these equations define the functions Fx and 
Fy up to arbitrary additive constants. 

The term with the square of the current j in Eq. (|20| ) 
can be split by the Hubbard-Stratonovich transforma- 
tion. 



2 



ihy 



2r' 



where b is the new auxiliary field defined at the links of 
the lattice. Although this field couples with the current 
like an additional gauge field, the integration in b is per- 
formed over all possible values including different gauges; 
this field is not a gauge field. By construction, the field 
b is defined only within the period of the system and, 
therefore, it satisfies periodic boundary conditions and 
does not introduce any non-trivial total magnetic flux. 
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Now let us resolve the conservation law (^l|) by intro- 
ducing an integer-valued scalar field h associated with 
the vertices of the dual lattice j^ = e^^'^V^h = V^/i. This 
definition implies that the increments of the additional 
field h across the system are related to the total currents 
in perpendicular directions, 



/h.M 



-I„ and 



A,/. 



I.\ 



(24) 



because of the periodicity and the current conservation 
these increments remain constant along the edges. 
At this point the effective action can be written as 



iSy 



b2 

— + i(b + A)V/i 



dS + iSsurf: 



(25) 



where the surface patching term Ssmi is given by Eq. (22) 



To arrive at the dual representation of the model we 
need to ^integrate the gauge coupling by parts, using the 
identitja 

A^e''''V^/i = AxV/i = /iVxA- Vx (A/i). 

The last term of this expression results in a contour in- 
tegral along the edge of the system, which can be split 
into vertical and horizontal parts, 

hAd\= I Ay{A^h)dx- I \^{Ayh)dy. (26) 
Jg Jo 

With the increments of fields A and h known, we can 
transform the increments of their products using the ex- 
pression 

A,{fg) = /(L„ y)A,5 + g{L,,y)AJ - A,/A,5, (27) 

so that the contour integral (|2l) takes the form 



hAdl^ / {V^Fyh + hA^~I^V^Fy)dx (28) 



(VyF^ h - lyAy + lyVyF^) dy . 



The first term in each of these integrals can be integrated 
by parts, with the integrated part once again evaluated 
using Eq. (|27|). With the expression for the total flux of 
the magnetic field, $ — A F^ — A.^Fy, following imme- 
diately from the Stokes formula and the definition (23), 
we obtain 



hAd\= / ixF^dy 




iyFydx - h(Lx,'Ly)^ 



A^dx + $-F^(Lj,) 



AydX — $ — FyiJj^ 



(29) 



*The product of the lattice fields Ah is defined by taking 
the field h at the plaquette shifted in the positive direction 
(up or right) of the link with the field A 



The first two terms cancel exactly with the edge patching 
term (p2|), while the remaining terms can be eliminated 
by properly defining the so far arbitrary additive con- 
stants for the scalar field h and the patching functions 
Fi. Specifically, defining /i(Lj;,Lj,) = 0, 



and 



F.Ci^y) 



Fy{Lx) 



Ax{x,'Ly)dx + $ 



Ay(L^,y)dy-^, 



eliminates the last two terms in Eq. (p9|), so that the 
performed integration by parts can be written simply as 



AxVhdS + 5surf = hS/ X Ads 



Although the described procedure works generally for 
any gauge field A, we had to be careful to separate the 
contribution from the Hubbard-Stratonovich field b: the 
gauge-fixing terms cannot depend on this field. The field 
b is periodic and the equivalent series of transformations 
is much simpler, immediately leading to 



hhdl^ L, 



hrdx + L, 



hydy 



(30) 



instead of Eq. ( p9| ) . It is convenient to use the potential 
representation b^ = Vif + e^^Vjg with the condition that 
the normal derivatives of both potentials vanish at the 
boundary, dg/dn = df /dn = 0. Then it follows from 
the periodicity 



h^dS 



{Vff + {Vgfdxdy, 



and the transformed action ( |25| ) can be written as 

^Sv=j[^+^h{Wy.A-y^g)yS+j^^dS 
+i{h + ly) /(L„ Lj,) - z/, /(O, Ly) - ily /(L„ 0), 

with the free terms in the last line resulting from the 
line integrals (BQ). It is clear that the two potentials / 
and g in this action are completely independent of each 
other. The summation over integer h only fixes the vortex 
charges v = V x A — V^y to integer values modulo 27r. 
Since the field b is single-valued, it does not contribute 
any additional flux and the total vorticity must be exactly 
equal to the total flux of the external magnetic field. 

Furthermore, the term J{Wg)^dS/2T precisely equals 
to the Coulomb energy associated with the vortex charges 
in the uniform background. This can be proven by intro- 
ducing an auxiliary field (f), 



(V5)2 T(V0)2 



2T 



\- iW(/)Wg 
T{V4>f 



Tiv^y 



i4N^g 



2-Ki(f)v + i(j)\7xA, 
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where the integrated part vanished because of the bound- 
ary condition dg/dn — 0. 

The considered transverse part of the dual action is 
independent of the lateral currents Ix , ly encoded in the 
discontinuity of the boundary condition (pj) for the field 
h. These currents act as the only sources of the scalar po- 
tential /. The energy of the equivalent Coulomb problem, 
quadratic in these currents, provides the only current- 
dependent contribution to the free energy; this contribu- 
tion is obviously non-critical. 



II. DUALITY IN THREE DIMENSIONS 

Now let us resume the analysis of the Chern-Simons 
gauge model in three dimensions. The major steps in- 
clude: (i) Start with the canonical form (|l|) of the model 
and rewrite it in terms of conserved integer-valued cur- 
rent j and auxiliary field b. (ii) Resolve the current con- 
servation law by introducing an additional integer- valued 
field m, j = Vxm, integrate the current coupling by parts 
and prove that, as a result of summation in m, the flux 
of the total gauge field through every plaquette is an in- 
teger, identified with the dual conserved current. Up to 
this point all transformations will be done exactly for any 
form of the original x-y model, (iii) In the Villain approx- 
imation, integrate over the longitudinal part of the field 
b to allow treating of its remaining part as a gauge field. 
(iv) Finally, resolve the constraint on the original Chern- 
Simons gauge field by yet another Hubbard-Stratonovich 
transformation, and express the result in the canonical 
form (0). 

i. Introduce the integer-valued current: Before intro- 
ducing the surface patching terms for the model (|l|) , it is 
convenient to transform its scalar part (g) to the equiv- 
alent current representation 



iSe = ^X^J ■ (Vfii-a-A-b) - ^ 



cosb 



(31) 



where |Lhe original form (0) is restored after the sum- 
mationH over all integer j and subsequent integration 
over the auxiliary field b at every link. As in the two- 
dimensional example, the integration over the phases 9 
imposes the conservation constraint V • j = 0, so that j 
can be identified as the integer- valued current. 



^The summation in discrete values of j automatically ac- 
counts for the periodicity of the cosine, validating its expan- 
sion in powers of b. For example, the quadratic approxi- 
mation cosb -^ 1 — b'^/2 is equivalent to the naive Villain 
approximation for the original model (H). The regular ex- 
pansion in powers of b can be constructed by analytically 
continuing the Fourier coefficients Im{T~^) of exp[(cos6)/r] 
to all real values of m, and using the function x{b) — 
Tin J JA(r~^) cos(A6)dA/(27r) instead of cosb in the expo- 
nent of Eq. (|l|). 



Assuming the kernel JC is transverse {i.e. (JCa) gauge- 
invariant,) it is easy to find the required form 



^^surf — ^ 



F, j 



-/Ca ) • ds 



(32) 



of the surface terms required to fix the gauge-invariance 
in the finite system. Here Fi are the patching functions 
for the sum of two gauge fields atot = a + A as defined in 
Appendix A, and the integration areas Si are the three 
visible faces of the cube as shown in FIG. |[ Since the 
auxiliary field b is single- valued, it creates no total mag- 
netic flux through the system and there are no boundary 
patching terms associated with this field. 



t^ c. 



X 



^^i^l 


^^tl 


1 
1 
1 
1 


.^ 



y 



FIG. 3. The integration areas Si in the surface patching 
terms (p3) are limited by the contours d, i — x, y, z. 



The surface patching terms ( |32| ) make the model 
gauge-invariant, but they also eliminate the integration 
over the constant part of the fluctuating gauge fleld a 
along with the resulting mean-field constraint 



ICadn — 



jdfl. 



(33) 



For the model with purely Chern-Simons gauge ker- 
nel (|g), Eq. (33) is the condition on the average flux of 
the associated magnetic fleld (V x a) = 27ra(j). This 
constraint allows us to associate the model (0) with the 
continuum CSLG model (^ and ultimately with the 
real quantum Hall systems; therefore the constraint (p3) 
needs to be imposed on the original model (|l|) by inte- 
gration over an additional uniform auxiliary field. 

ii. Resolve the current conservation: By analogy with 
the two-dimensional problem, we can resolve the lattice 
current-conservation constraint V • j = by introducing 
an auxiliary field m. 



j = Vx m. 



(34) 



There is, however, an important difference with the two- 
dimensional case: in 3D the gauge freedom associated 
with the solution of Eq. (|3J) is much more broad. Simul- 
taneously, in generic gauge m is not an integer-valued 
field, and even the existence of the integer-valued solu- 
tion is not immediately obvious. 

Let us choose the gauge m^ = 0, then it is easy to 
check that the sums 
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■m^{x,y,z)= I jy{x,y,C)dC~ / jzity,0)dC, 
Jo Jo 



my{x,y,z) 



Jx{x,yX)dC, 



(35) 



satisfy Eq. (pj); by construction this solution is integer- 
valued. Now we may introduce the patching functions 
Gi associated with the integer-valued gauge field m. Ex- 
plicit calculation shows that both Gx = Q and Gy = up 
to an integer constant, and 

dxGz{x,y)^ / jydz, dyG^ix.y) ^ ~ / jydz\ 

Jo Jo 

these equations have a solution since the curl 



dyid^Gz) - dxidyGz) = / {d^jx + dyjy) dz 
Jo 



dzjz dz = 



vanishes. According to Eq. (Al), the complete current- 
dependent coupling can be integrated by parts as 

/ V X m • (b + atot) dn + Ssuri = / m • V x (b + atot) dn 
+ J2 hf^h,dx,+jG,[\7x(h + atot)lds, (36) 

i—x,y,z ^ 

where atot = a + A. Independent summations over rux 
and niy restrict the x- and y- components of the field 



v= --V x(b + atot) 

ZTT 



(37) 



to integer values. This field is a total curl, and conse- 
quently it is conserved, V-v = 0. Therefore, the dif- 
ferences dz^^z of the z-component of this field at adjoin- 
ing links are also integer-valued. Finally, the summation 
of the patching surface term at the surface Sz over the 
values of the function Gz restricts the values of v^ at 
this surface to integer values, completing the proof that 
all components of the field v can take only integer val- 
ues. This conserved field v will be identified as the dual 
integer- valued current. 

Up to this point all transformations were exact and 
very general, valid for virtually any imaginable form of 
the lattice gauge model. Eq. ( |3^ ) relates the lines of the 
dual current to the flux lines of the combined magnetic 
field Vx(A-f a-f b), so that the total dual current through 
the system 



1 /■ * 

V ds = — / V x (A + a) ds = — (38) 

2tt I 2tt 



in any given direction is precisely equal to the total num- 
ber of flux quanta $/27r of the combined magnetic field 
in this direction (the contribution of the periodic auxil- 
iary field b vanishes after the surface integration). This 



universal relationship between the total current and the 
total magnetic flux in mutually dual models is one of the 
central results of this work. Previously, a similar result, 
relating the flux of the magnetic field witlt-tlie number of 
vortices, was obtained by Lee and Fishertili23, in an ap- 
proximate lattice duality involving a "softening" of the 
integer constraint on the total number of vortices. 

iii. Integrate away the scalar potential: The goal of 
the subsequent steps is to rewrite the dual action in a 
form possibly close to that of the original action (g). This 
can be done only approximately, ignoring some terms ir- 
relevant in the RG sense. In particular, the integration 
over the auxiliary gauge field b will be done only in the 
Villain approximation, leaving out the non-linear correc- 
tions. With this in mind, the result of the summation 
of the combined action (0) , (p^) over the integer- valued 
current j can be written as 

/■ b^ 1 /■ 

J 2T 2 J 

+i 2_] ^p- / bp dx^ + (surface terms) . (39) 

fj.—x,y,z ^ 

Here the field a is a solution of the equation 

Vxa=27rv- Vx(b-FA), (40) 

while the Hubbard-Stratonovich field b and the con- 
served dual integer- valued current v with given total flux 
J -V ds — $/27r through the sample are independent vari- 
ables. Thepresence of surface terms ensure that the ex- 
pression (|3^) does not depend on the particular gauge 
chosen for the solution a of Eq. (|4(]|). 

Following the two-dimensional example, we can elimi- 
nate the integrals of b in the second line of Eq. (p^ ) by 
writing b = Wcj) -|- V x g = V(^ -I- b', with the tangential 
part of the second term vanishing everywhere at the sur- 
face of the system. Then the integration over the field 
(j) yields a non-critical quadratic in the currents I^ con- 
tribution to the free energy, proportional to the energy 
of the Coulomb interaction of the charges ±1^ set in the 
corners of the cube in FIG. and their mirror images. 

Now the constraint on the boundary values of the 
field b' can be removed by adding the appropriate sur- 
face patching terms rendering the expression gauge- 
independent; this is possible because the surface terms 
evaluate to zero with the original boundary condition 
of vanishing tangential components of the field b'. The 
gauge fields a and b' enter Eq. (B3) only in the combi- 
nation a' = a -f b'; integrating awav-,their difference we 
obtain the combined gauge couplingCJ 

M! + ia/Ca^ ia'/C'a', /C' = {T + IC-^' , 

where the field a' satisfies the equation 

Vxa' = 27rv- VxA, 

and all inverse matrices must be evaluated in the trans- 
verse gauge V • a = 0. 
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iv. Introduce new auxiliary field: Finally, we need to 
express the coupling a'/C'a' as a path integral over an 
additional auxiliary field c. With the appropriate sur- 
face patching integrals implicit, this can be done by the 
Hubbard-Stratonovich transformation 

ia'/C'a' ^ ia' • V X c + i (V X c)/C'"^ (V x c), (41) 

where the additional constraint on the total flux of the 
auxiliary field through the system 



V X c ds = — i / IC'sL ds = —I 



(42) 



is mandated by the degeneracy of the gauge kernel JC' at 
zero momentum due to the presence of the surface patch- 
ing terms. Indeed, the Hubbard-Stratonovich transfor- 
mation (|4l| ) can be interpreted as the addition in the 
exponent 

a'/C'a' ^ a'/C'a' + (V x c + i/C'a') /C'"^ (V x c + i/C'a') , 

with the subsequent integration in a'; the condition (E2) 
is required for the inverse operator to be meaningful. 

With the help of Eq. (Eoh, the total Lagrangian be- 
comes 

;r , .^ (Vxc)/C'"VVxc) , , 

After defining the dual gauge field 

V X a = 27rV x c - 27ri JC'A, 
we can write the Lagrangian of the system in the form 



iv(a-27r(iVx)"'/C'A) 



where the dual gauge kernel 

P x /C'"^ X P rp2 p X /C-i X P 



(44) 



JC 



(2^) 



f27r) 



(2^)2 



and the total flux $ of the gauge fields coupled to the dual 
current v is determined by the original particle current. 



$ = / (V X a + 27ri/C'A) ds = -2tiI 



(45) 



as immediately follows from Eq. ( p^ ) and the definition 
of the dual gauge field a. It is easy to see that the second 
duality transformation restores the form of the model up 
to the trivial symmetry (|^). 

In the limit of small momenta or slow- varying external 
field A, and assuming the original gauge kernel (^) has 
purely Chern-Simons form (0), the important part of the 
dual model can be written as 



£(^^^) 



*i^(Vx&)^ (46) 



indicating that, up to irrelevant difference in the quasi- 
particle's size, the dual model is the Chern-Simons Vil- 
lain x-y model with the excitations of electric charge 
1/a, statisticsL — 1/a, and the dimensionful lattice 
temperatureCjO T = Tjc? . 

The global conditions on the total dual fiux (|45|) and 
current ( p8|) restore the mean-field constraint (|33|) for the 
dual model only if the gauge kernel of the original model 
has purely Chern-Simons form (0). Indeed, 



Vxads 



© 



* 



27r 
a 



a 



-27rl- 
27r 



Ol 



1 



* 



a 



where ^a denotes the total fiux of the external magnetic 
field. It is this property that shows that the chosen def- 
inition of the gauge coupling in finite systems is natural 
and corresponds to the expected behavior for the model 
describing quantum Hall systems. The relationship of 
the considered class of models to the quantum Hall effect 



will be discussed more specifically in Sec. IV 



Unfortunately, it is not clear how to perform the transi- 
tion between Eqns. (E^) and (E3) for the generalized farm 
of the Chern-Simons coupling (^, which is supposedEZi to 
describe the even-denominator (metallic) fixed points of 
the quantum Hall system. In this case the background 
charge density fluctuations in Eq. (E^) cannot be self- 
consistently traded for the modified magnetic field in 
Eq. ( [44| ) , indicating that additional entities like the "mag- 
netic" chargesE3 are needed. At this time, it is not clear 
what physical meaning can be attributed to such objects. 



III. FLUX ATTACHMENT TRANSFORMATION 

In addition to the duality, the continuum Chern- 
Simons models have another important symmetry, the 
periodicity, or the symmetry under the fiux attachment 
transformation. In this section I show that the lattice 
version of this transformation is exact aad, in addition 
to modifying the Chern-Simons cocfficicntcj, changes the 
magnetic flux $ by the amount proportional to the global 
current I. 

Let us fix a configuration of integer- valued conserved 
current j and define the related gauge field Cj as a solution 
of the equation V x Cj = 47rj. In a system where the 
current lines never cross the boundary, the coupling 



X 



Cj j dil + surface terms 



(47) 



is always an integer times 27r. Because of the surface 
patching terms, this coupling is always an integer times 
27r, independent of the gauge chosen for the field Cj. In- 
deed, for the current lines that do cross the boundary, 
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the surface terms can be interpreted as the "back flow" 
currents, as shown in FIG. |l], and the statement for the 
closed current lines apply. 

We would like to rewrite the coupling (^Tj) as a path 
integral over a dynamical variable c instead, 



X = 



cj 



(c9c) 
87r 



dri + surface terms. 



(48) 



but the explicit calculation shows that this expression 
does not contain an integration variable needed to fix 
the total flux of the fleld c. To set this flux correctly, we 
need to add the term 



A„ 



Vxc 
47r 



with integration over the components of the uniform aux- 
iliary fleld A^. This is equivalent to restoring the inte- 
gration over the uniform part of the fleld c removed by 
the surface patching terms; the whole expression still re- 
mains gauge- invariant. With the addition of this con- 
straint, the integration in Eq. (Bsf) yields exactly the self- 
couphng (p^l). 

Because it is always an integer multiple of 2tt, the ex- 
pression ( p8| ) can be introduced as an additional phase 
to the partition function ffl) without changing its value. 
As a result, the full partition function contains a mini- 
mal coupling of the integer- valued current j to the sum of 
two fields a' = a -|- c with their respective gauge kernek 
/C^'^ and K.^" = ie'^^PV p/{Att). After integrating awayN 
the difference of these fields, the resulting model has the 
gauge-invariant coupling to the fleld a' with the gauge 
kernel 



/C'= [IC-^+JC^^] 



(49) 



with the Chern-Simons coefficient a' = a -I- 2. Since 
both fields a and c were constrained to have a fixed total 
magnetic flux, the combined fleld a' has the magnetic 
flux 



$' = * + 47rl 



(50) 



in every direction. Clearly, the total magnetic flux in- 
creased by twice the number of the particle's world lines 
propagating through the system in the same direction, 
while the total particle current I remains unchanged. 

Similarly to the duality, the flux attachment trans- 
formation preserves the mean-field constraint (p3) for 
the models with purely Chern-Simons gauge kernel (ra). 
Therefore, imposing such a constraint once for the initial 
model immediately sets it for the whole hierarchy of mod- 
els generated by the duality and flux attachment trans- 
formations, in agreement with our expectations for the 
models describing plateau transitions in quantum Hall 
systems. 



IV. CORRESPONDENCE LAWS 

The hierarchy transformations of the CS x-y model 
change the microscopic dimensionless charge a deflned 
through the gauge coupling kernel (|4|), and, more im- 
portantly, they also modify the global parameters, the 
Euclidean current and the magnetic flux, in a univer- 
sal fashion. It turns out that these global transfor- 
mations are formally equivalent to the Correspaadence 
Laws, suggested^ and conflrmed experimentallyBESa'I'Q, 
as the symmetries of the quantum Hall phase diagram. 
This analogy can be seen in the real-time representation, 
which requires treating the charge-density p and the lat- 
eral current j separately. 

i. Symmetries of the phase diagram. The Wick rota- 
tion transforms the third component jz of the Euclidean 
current into the charge density p, while the corresponding 
component of the total magnetic field Btot — Vx (a+ A) 
is transformed to the component of the magnetic field 
normal to the two-dimensional plane of motion of par- 
ticles. It is convenient to express the global transfor- 
mations of the x-y model in terms of the bosonic filling 
fraction v^ = 2tt{p) / (Btot)- 



reflection (|| 

duality (H), (@ 

flux attachment (^ 



.,-1 



-t^, 



(51) 
(52) 
(53) 



which, are formally equivalent to the Correspondence 
Lawaa for fermions. Indeed, by attaching a flux quan- 
tum to every boson, we can convert it to a fermion in 
a modifled magnetic fleld, their flUing fractions related 
through the identity l/i^ = l/t'f — 1. Then each of the 
basic transformations for fermions become equivalent to 
a combination of transformations (pl|)-(p3[), as explicitly 
shown in Table H. 



Transformation 



Landau level addition 
flux attachment 
particle-hole symmetry 




Bosons 



f b 



fb 



2 
2 



-b-^ 



TABLE 1. The Laws of Correspondence in terms of Fermi- 
onic ut and Bosonic Vh filling fractions related by the identity 

1/Ub = IjVi - 1. 

ExperimentallyoHu, the duality symmetry was ob- 
served for the transitions between the quantum Hall liq- 
uid and insulator. Speciflcally, the 2DEG demonstrates 
reciprocal longitudinal resistance and the mirror symme- 
try of non-linear I-V curves in the points at the phase 
diagram, related by the expression 



{a-v-^){c 



z>-M 



= 1, 



(54) 



where a is an odd integer value of the quantized Hall 
resistance in the corresponding plateaus, measured in 
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the units of h/e^, and v = 2-Kp/B is the measured fill- 
ing fraction of the electrons, which is varied across the 
transition. In terms of the "bosonic" filling fractions 
i/b = (a — v^^)^^ , this relationship can be written sim- 
ply as z/b I'b = 1, which is equivalent to the combination 
of the duality (|5^) with the reflection (pl|). Note, that 
the measured value of the filling fraction Vc in the transi- 
tion point is close but not exactly equal to the theoretical 
value V* = (a-f 1)^^, implying that the quasiparticles are 
not in one-to-one correspondence with the original elec- 
trons. A partial mixing of the opposite spin states, not 
included in the theoretical model, provides a qualitative 
explanation for this discrepancy. 

ii. Symmetries of transport coefficients. The remain- 
ing two components of the Euclidean current j = (p, if) 
originate from the lateral current J of the 2+1 dimensional 
system, while the in-plane Euclidean magnetic field is 
generated by the electric field 



B, 



iE; 



v 



B, 



-lE:r^ 



(55) 



The rotation by 7r/2 is needed to transform the com- 
ponents of the axial vector B into the polar vector E; 
the inversion of the handedness changes both the sign of 
B and the sign of the rotation angle, so that the rela- 
tion ( |55| ) remains valid. Unlike that for the usual charge 
density and magnetic field, the Wick rotation of spatial 
components involve the complex phase i = e*'^/^, indicat- 
ing that the analytic continuation in the values of these 
fields is needed. Any singularity encountered during this 
procedure would invalidate the current -voltage mapping. 
This is certainly related to the fact that the presence of 
the electric field E and the current density J may involve 
dissipation. However, this mapping is valid in all orders 
of the perturbation expansion in powers of the electric 
field E; the actual region of validity of this procedure is 
determined by the convergence of the perturbation se- 
ries, as well as the presence and the relative importance 
of non-perturbative terms. 

The duality interchanges the components of the cur- 
rent density with the appropriate components of the ro- 
tated electric field. In terms of physical electric field and 
current density, the duality relationship can be written 
as 



Er, 



2^5y,: 



±a 






±« -1, 



(56) 



where the first (second) indices correspond to top (bot- 
tom) signs. In the linear regime the substitution E = pj 

and E = p] shows that the equation (|56| ) is equivalent to 
the inversion of the bosonic resistivity tensor 



Pb = (^b 



o-fe = Pb 



(57) 



between the mutually dual points. Here the bosonic resis- 
tivity tensors are related to the usual resistivity tensors. 



measured in the units of hc/e — 2'Kfic/e, by the relation- 
ship 



Pb = p- a 



1 
-1 



(58) 



As pointed out in Rcf. E3, Eq. ( pT\ j corresponds to the 
inversion of the longitudinal resistivity p^x = ^/ Pxx ob- 
served in experimenta as long as the Hall resistivity re- 
mains equal to its quantized value pxy = a, or, equiva- 
lently. 



2 I 2 
xy ' XX 



a aj. 



This statement, known as the Semicircle Law, is associ- 
ated with the presence of disorder, . and . is known to hold 
exactly for two-component systema23i2il. 

| iii . .Mniversality of phase transitions. Perturbative- 
lylljll3'Ea, the universality of phase transitions in Chern- 
Simons models does not seem to be a wsobable outcome. 
However, the universality was provencJ for the transi- 
tions in the clean CS x-y model where the "relativistic" 
symmetry is not broken by external magnetic field or 
total charge density. In the absence of these fields the 
independent parameter v is absent, and the hierarchy 
mapping between the vicinities of phase transitions at 
different values of a is complete. The exact universal- 
ity class depends on whether the given fractional value 
of the Chern-Simons coefficient originates from bosons 
or fermions, but it is not likely to correspond to that of 
any phase transition experimentally observed in quantum 
Hall systems. 

In the extended CS x-y model (P, where both the 
external magnetic field and the charge density can be 
present, there are at least two independent parameters 
that characterize the state of the equilibrium system: 
the filling fraction ly and the Chern-Simons coefficient 
a. Transforming both of them at each step, it is im- 
possible to construct a mapping between two quantum 
Hall states of electtpns {a = 1). Even though the proof 
of the universalitycj remains formally valid, in the pres- 
ence of non-zero filling fraction the universality applies 
not to the original electrons, but to vortices, particles 
with the transformed Chern-Simons coefficient, at the 
dual filling fraction. In the original formulation of the 
hierarchy theoryd, this problem is resolved by making 
an assumption about the independence of the bosonic 
current-current correlation function of the Chern-Simons 
coefficient. So far, even though an ample experimental 
evidence confirming this non-trivial assumption is avail- 
able, a rigorous theoretical proof, or even an effective 
model that would display this behavior, is missing. 

Any possible scenario explaining this puzzle must 
include the disorder. For example, if the disorder- 
associated scattering turns out strong enough to screen 
the correlations associated with the Chern-Simons inter- 
action at large distances, the effective disorder-averaged 
model will be independent of the Chern-Simons coeffi- 
cient. Since both the duality and the fiux attachment 
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preserve the shape of the disorder potential, it is con- 
ceivable that the disorder-averaged model would never- 
theless inherit the original hierarchy symmetries, at least 
in some range of parameters. This poir 
ported by the bulk of experimental dataMlii 
which indicates that the fractional and integer quantum 
Hall transitions are in the same universality class, and 
their vicinities are at least approximately related by the 
Correspondence Laws. 



ssasii 



In this context, a similar and|-.sn far unexplained 
current-voltage symmetry observedEfQ in the conductor- 
insulator transition in silicon MOSFETS at zero mag- 
netic field appears especially suggesting. Certainly, there 
remains a more mundane possibility that the observed 
symmetry is a mere consequence of the universal critical 
scaling. 
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CONCLUSION 

The Chern-Simons x-y model with external fields con- 
sidered here was constructed as a lattice regularization of 
the CSLG model. Even though this model is written in 
rotationally (relativistically) covariant form (|l])-(||), the 
non-zero charge-current density and external magnetic 
and electric fields break the Lorentz symmetry in exactly 
the same way as it is broken in the original continuum 
model. In addition, the studied model demonstrates the 
same non-local symmetries as the CSLG model — only at 
the lattice these symmetries can be formulated as exact 
statements, without the limitations of the RPA approxi- 
mation used in the continuum model. 

Unfortunately, neither model clears the issue of the 
observed universality of phase transition. Note, however, 
that both the duality and the flux attachment transfor- 
mation remain exact in the Chern-Simons x-y model even 
in the presence of disorder; the latter is merely carried 
from one representation to another without any signif- 
icant -change of the statistical properties. It is easy to 
checktJ that in the replicated version of this model the 
disorder averaging commutes with the duality transfor- 
mation. This gives a hope that there exists an effec- 
tive description of diffusion in this strongly-interacting 
system, independent of the Chern-Simons coupling and, 
therefore, symmetric with respect to hierarchy transfor- 
mations. 

Apart from this issue of universality, the analyzed CS 
x-y model appears to be a perfect candidate for provid- 
ing an effective description of the quantum Hall phase 
transitions. In this second-quantized model the non-local 
duality and flux attachment transformations can be per- 
formed exactly. These transformations are expressed in 
terms of system-averaged charge and current density and 
also the electric and magnetic fields; they generate the ex- 
perimentally observed global symmetries of the phase di- 
agram known as the Correspondence LawaJ, and the non- 
linear pcurrent-voltage mapping between mutually dual 
pointsEI. The global nature of the transformed parame- 
ters imply that they are not subject to renormalization. 
Moreover, the derived duality relationship is so robust 
and insensitive to particular details of both interaction 
and disorder, that the immediately following I-V reflec- 
tion symmetry could be viewed as a signature of the du- 
ality of this kind. 
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APPENDIX A: HERMITICITY OF THE 
CHERN-SIMONS COUPLING 

Here I derive the conditions of Hermiticity 

{adb)dn -J2[ ^''^ X b • ds (Al) 

= / {bda)dn - XI / -PiV X a • ds 

for the gauge-invariant Chern-Simons coupling in finite 
systems, with the patching functions Gi and Fi defined 
(up to the additive constants specified below) by the ex- 
pressions 

VGi = A,a, VFi = A,b (A2) 

Consider the difference of two volume integrals 

J= /aVxbdO- /bVxadO 

~ V ■ {h X a.) dil = (p h X a- ds 

transformed to the integral over the full outer surface 
of the system, which can be split into the contributions 
from three pairs of parallel faces, 

i={x,y,z} 



J,^ I VF, X (a--VGj -ds 

I VG, X (h-lvF^ -ds. 



(A3) 



In each term the integration is limited to the front Sx, 
right Sy or top S^ sides of the system as shown in FIG. ||, 
and the patching functions ( A2) were used to simplify the 
increments of the gauge fields across the system. Each of 
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these six integrals can be further transformed using the 
vector identity 

Vi^xa = Vx(aF)-FVxa 

into the sum of the contour integral and the surface in- 
tegral 



Some additional simplification can be achieved by using 
the expression for the total flux of the gauge field 



A^F. 



A,F, = 



a • dl = *e 



c^ 



VG' 



V_G 
2 



in z-direction, and similar expressions for other compo- 
nents of a and the field b. After some rather tedious but 



VF X a — \ ds — f F [a. — 1 • dl — / FVxa- ds. straightforward algebra the expression ( |A6| ) simplifies to 



In every case the surface integral cancels with one of 
the surface patching terms in J. Therefore, the Chern- 
Simons coupling with appropriate edge patching terms is 
a Hermitian operator as specified by Eq. (Al) if and only 
if the sum of all line integrals 



L 



E 



i^. (a-ic, 



dl 



a '^ b 



(A4) 



vanishes. 

The closed line integrals in Eq. (A4) can be simpli- 
fied by grouping together the contributions of the parallel 
parts of each contour, and once again using the patching 
functions for the differences of the gauge fields 



L 



a^y^ 



dy<F, 



(Fz - KFz 



:G, 



^G 



Gx,v + -2^x^z,y 



a ■ 
G- 



b 
F 



(A5) 



Here we imply the summation over all permutations of 
the coordinates, Gr 



dyGz, and the functions are eval- 



uated at the segment (La,,!/, L^) or its equivalents gen- 
erated by cyclic permutations of the coordinates. It i s 
convenient to use the symmetries of the expression ( [Aq ) 
and write symbolically 



L = 



12 



-A F, 



(F, - A,F,) G,^ 



:G, 



-A.^Gz^y 



Gx,y >dy 



(A6) 



implying the summation over six permutations of the co- 
ordinates and two additional permutations of the gauge 
fields a and b with their corresponding parity factors. 
This way one does not need to write all terms explicitly 
in order to simplify the expression (A6). For example, 
the integral 

J2 [ P-G-^y '^V^\Y. I (P-'^-.y + GxFz^y) dy 

A,(F,G,) 
^IJ2 A,yF, G, + Fz A^G, - AyFz A^G,. 



12 



12 

i5:(F.G.)i^.i^, 

12 12 



F-T.^t(Fy^ 



bydy 



.^b 



.^b 



(A7) 



The function Fy in this expression is evaluated at the 
vertex (L^,, Lj,, L^) and the integration in y is performed 
along the line a; = L^, z = L,,. 



With the help of Eq. ( A7) we can fix the so far arbitrary 
constants in the definition ( |A2| ) of the surface patching 
functions Fi and G^. Specifically, by defining 

Ly 

Fy{L,,Ly,U) ^ lby{L^,y,Lz)dy^^^l + i*^ (A8) 



(and similar definitions for other functions Fi and Gi) 
we can completely suppress the contour integrals L and 
therefore render the Chern-Simons coupling with the sur- 
face term a Hermitian operator as defin ed by Eq. (Al). 
It is important that the condition (A8) is not modified 
by the gauge transformation (n^ because the fluxes ^^ 
are gauge- invariant. Therefore, both the Hermiticity and 
gauge-invariance can hold at the same time. 
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